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Abstract 

We investigate a steady flow of compressible fluid with inflow boundary condition on the 
density and slip boundary conditions on the velocity in a square domain Q E R 2 . We show 
existence if a solution (v, p) E Wp (Q) x W p l (Q) that is a small perturbation of a constant flow 
(v = [1,0], p = 1). We also show that this solution is unique in a class of small perturbations 
of the constant flow (v, p). In order show the existence of the solution we adapt the techniques 
know from the theory of weak solutions. We apply the method of elliptic regularization and a 
fixed point argument. 
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1 Introduction and main results 



The problems of steady compressible flows described by the Navier-Stokes equations are usu- 
ally considered with the homogeneous Dirichlet boundary conditions on the velocity. It is worth 
from the mathematical point of view, as well as in the eye of applications, to investigate different 
types of boundary conditions. A significant feature of the compressible Navier-Stokes system 
is its mixed character: the continuity equation is elliptic in the velocity whereas the continuity 
equation is hyperbolic in the density. If we assume that the flow enters the domain, then the 
hyperbolicity of the continuity equation makes it necessary to prescribe the density on the inflow 
part of the boundary. A time-dependent compressible flow with inflow boundary condition has 
been considered by Valli and Zajaczkowski in [|2T|. The authors showed existence of a global 
in time solutions under some smallness assumptions on the data. They also obtained a stability 
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result and existence of a stationary solution. Plotnikov and Sokolovski investigated shape opti- 
mization problems with inflow boundary condition in 2D [16J and 3D |[T5ll . working with weak 
solutions. Regular solutions to problems with inflow boundary conditions have been investigated 
mainly by Kweon in a joint work with Kellogg and with Song BH. The results obtained by 
these authors require some assumptions on the geometry of the boundary in the neighbourhood 
of the points where the inflow and outflow parts of the boundary meet. In [[8]| Kweon and Kellogg 
investigated the case when the inflow and outflow parts of the boundary are separated, obtaining 
regular solutions. What seems to be interesting is to investigate an inflow condition on the den- 
sity combined with slip boundary conditions on the velocity, that allow to describe precisely the 
action between the fluid and the boundary. The slip boundary conditions have been investigated 
by Mucha iTTOll for incompressible flows, and also by Fujita [4] and Mucha and Pokorny ifTTI for 
compressible flows. 

Here we investigate a steady flow of a viscous, barotropic, compressible fluid in a square 
domain in R 2 satisfying inhomogeneous slip boundary conditions on the velocity combined with 
an inflow condition on the density. We impose that there is no flux across the bottom and the 
top of the square, so that it can be considered a finite, two dimensional pipe. From the analytical 
point of view our domain prevents the singularity that appears in a general domain where the 
inflow and outflow parts of the boundary coincide. 

We show existence of a solution that can be considered as a perturbation of a constant solution 
(v = (1,0), p = 0). Under some smallness assumptions we can show an a priori estimate in 
a space Wp(Q) x Wp(Q) that is crucial in the proof of existence of the solution. Now let us 
formulate the problem under consideration more precisely. 

The stationary compressible Navier-Stokes system describing the motion of the fluid, sup- 
plied with the slip boundary conditions, reads 

pv ■ Vf — pAv — (p + v) Vdiv v + Vp{p) = in Q, 

div (pv) = in Q, 

n-T(v,p^)-r + fv-r = b on T, (1.1) 

n ■ v = d on T, 

P = Pin on r in , 

where Q = [0, 1] x [0, 1] is a square domain in R 2 with the boundary T and T in = {x E T : 
v ■ n(x) < 0}. We will also denote T out = {x E T : v ■ n(x) > 0} and r = {x E T : 
v-n(x) = 0}. Next, b E Wp~ 1/p (T),dE Wp~ 1/p (T) and p in E Wp~ 1/p (T in ) are given functions. 
v = (v^\ v^) is the velocity field of the fluid and p is the density of the fluid. We assume that 
the pressure is a function of the density of the form p(p) = p 1 for some 7 > 1. The outward unit 
normal and tangent vectors are denoted respectively by n and r. We assume d = on T , what 
means that there is no flow across these parts of the boundary. Moreover, 



is the stress tensor and 



T(v, p) = 2pD(v) + v div v I — pi 
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is the deformation tensor. \i and v are viscosity constants satisfying \i > and v + 2p > 
and / > is a friction coefficient. The slip boundary conditions (ll.llh /i are supplied with the 
condition (11 .Ilk prescribing the values of the density on the inflow part of the boundary. Under 
the assumptions on /i and v the momentum equation (ll.llh is elliptic in u, whereas the continuity 
equation (11.11) 9 is hyperbolic in p. 

Our method would also work with no modification if we considered a perturbation of the 
constant flow (v, p) satisfying (ll.ll) i with a term pF on the r.h.s provided that | \F\ \ L was small 
enough. 

Since T(v, p 1 ) = 0, the constant flow (v, p) fulfills equations (11.11) with boundary conditions 
fv-r = b — fr^ and n ■ v = d — fr^\ 
Our main result is 

Theorem 1. Assume that \\b — fr^\\ w i-i/ P , rv \\d — n^'ll 2-i/ Pfr , and \\p in — 111 . are 

small enough and f is large enough. Then there exists a solution (v, p) G Wp (Q) x WhQ) to 
the system < \1.1\) and 

\\v-v\\w* + \\p-p\\w2<E, (1-2) 

where E is a constant depending on the data, i.e. on d, pi n , b, the constants in the equation and 
the domain, that can be arbitrarily small provided that the data is small enough. 

Moreover, if (vi, pi) and (v 2 , P2) are two solutions to M.l\) satisfying the estimate A1.2\) then 

There are several difficulties in the proof of Theorem [T] that result, roughly speaking, from the 
mixed character of the problem. In a general domain a singularity appears in the points where the 
inflow and outflow parts of the boundary meet and we can not apply the method used in this paper 
to obtain an a priori estimate. However, there is another difficulty in the analysis of the steady 
compressible Navier-Stokes system, independent on the domain. This difficulty lies in the term 
u ■ Vw. Namely, if we want to apply some fixed point method then this term makes it impossible 
to show the compactness of the solution operator. We overcome this difficulty applying the 
method of elliptic regularization. We solve a sequence of approximate elliptic problems and 
show that this sequence converges to the solution of (11.11) . This is a well-known method that has 
been usually applied to the issue of weak solutions (" lfT4l . IfTTlO . and differs from the approach of 
Kweon and Kellogg used to derive regular solutions in Q, jH. 

Let us now outline the strategy of the proof, and thus the structure of the paper. In section[2]we 
start with removing inhomogeneity from the boundary conditions (ll.llh a. It leads to the system 
(12.31 ), and we can focus on this system instead of (11.11) . In the same section we define an e - elliptic 
regularization to the system (12.31 ) and introduce its linearization (12.4b . In section|3]we derive an e - 
independent estimate on a solution of the linearized elliptic system (Theorem[2]). Although linear, 
the system (12.41) has variable coefficients and thus its solution is not straightforward. In order to 
solve (12.41) we apply the Leray-Schauder fixed point theorem in section HI using a modification 
of the estimate from Theorem |2l In section [5] we use the a priori estimate to apply the Schauder 
fixed point theorem to solve the approximate elliptic systems. In section [6] we prove our main 
result, Theorem \T\ The proof is divided into two steps. First we show that the sequence of 
approximate solutions converges to the solution of (12.31) and thus prove the existence of the 
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solution to (11.11) satisfying the estimate (11.21) . Next we show that this solution is unique in a class 
of small perturbations of the constant flow (v, p). We see that the estimate from Theorem[3]is in 
fact used at three stages of the proof, therefore we show it in a detailed way in section [3l 

2 Preliminaries 

In this section we remove the inhomogeneity from the boundary conditions (ll.lb /i 5. Then we 
define an e - elliptic regularization to the system (OQ). We also make some remarks concerning the 
notation. Let us construct u £ W^(Q) and w £ W^(Q) such that 



n ■ w |r 



d - n {1) and w \ Tin = p in - 1. (2.1) 



Due to the assumption of smallness of d — rv- 1 ' |r and pi n — l|r in we can assume that 

||«o||w2, ||wo||w| << 1- (2.2) 

Now we consider 

u = v — v — Uq and w = p — p — wo- 
One can easily verify that (u, w) satisfies the following system: 

d Xl u — pAu — [y + p)Vdivu + j(w + wq + l) 7 "^^ = F(u, w) in Q, 

(w + wo + 1) div u + d Xl w + {u + uo) ■ Vu> = G(u, w) in Q, 

n ■ 2pB(u) -t + f u-t = B onT, (2.3) 

n ■ u = on T, 

w = on T in , 

where 

F(u, w) = — (w + w + 1) (u ■ Vu + u ■ Vm ) — w (u ■ Vit ) 
— (w + wo + 1) u ■ Vw — 7(10 + wo + 1) 7_1 Vwo + pAu + {y + p)Vdiv u Q — (w + l)u ■ Vu , 
G(u, w) = — (w + wo + 1) divuo — {u + uo) ■ Vwo — d Xl wo 

and 

B = b - 2pn ■ B(u ) ■ r - /r (1) . 

In order to prove Theorem[T]it is enough to prove the existence of a solution (u, w) to the system 
(12.31) provided that ||uo||w|j H^ollvKp 1 and ||S|| w i-i/p, r ^ are small enough. As we already men- 
tioned, the presence of the term u ■ Vw in the continuity equation makes it impossible to show the 
compactness of a solution operator if we try to apply fixed point methods directly to the system 



(12.31) . We overcome this difficulty applying the method of elliptic regularization. The method 
consists of adding an elliptic term — eAw to the r.h.s of (12.3b -? and introducing an additional 
Neumann boundary condition. Since the density is already prescribed on the inflow part of the 
boundary by (|2.3I) r, we impose the Neumann condition only on the remaining part of the bound- 
ary. While we are passing to the limit with the density in Wp - norm, the Neumann condition will 
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disappear. Similar approach has been applied to the issue of inviscid limit for the incompressible 
Euler system in @. Consider a following linear system with variable coefficients: 



d xi u e — fiAu e — {y + ji)Vdivu e + j(w + wo + l) 7_1 Vw e = F e (u, w) in Q, 

(w + w + 1) div u t + d Xl w e + {u + Uq) ■ Vw e — eAw t = G € (u, w) in Q, 

n ■ 2/iD(w e ) • t + f u e ■ r = B on T, 

n ■ u e = on T, 



(2.4) 



ito = on r \ r 



171) 

in- 



where (u,w) G x are given functions and F e (u,w) and G e (u,w) are regular- 

izations to F(u, w) and to) obtained by replacing the functions u and w by their regular 
approximations u e and w e Q . 

Let us define an operator T e : V c W%{Q) x H^(Q) -> W p 2 (Q) x W*(Q) : 

(« e , w £ ) = ^(m, w) <^=^» (w e , w e ) is a solution to (12.41) . (2.5) 

where P is a subset of W% (Q) x that we will define later. Using the operator T e we 



define an e - elliptic regularization to the system (12.31) . 

Definition 1. By an e - elliptic regularization to the system ( 12.31) we mean a system 

(Me, We) = T e («e, We). (2.6) 

We want to show the existence of a solution to the e - elliptic regularization to the system (|2.3I) 
applying the Shauder fixed point theorem. The strategy has been outlined in the introduction. In 
section|4]we show that T e is well defined, which means that for given (u, w) there exists a unique 
solution to (12.41) (Theorem [3]). In fact we show that T t is well defined for e small enough, but it 
suffices since we are interested in small values of e. 

In section [5] we show that T e satisfies the assumptions of the Schauder fixed point theorem 
and thus we solve the system (12.61) for e small enough. 



As we already said, the key point is to derive an e - independent estimate for the system (12.41) . 
which is used at different stages of the proof. We derive such estimate in the next section. Before 
we proceed, we will finish this introductory part with a few remarks concerning notation. 

For simplicity we will denote 

7(w + w + l) 7 
«oW = — , 

ai (w) = 7 (w + w + l) 7 " 1 , 1 ; 

a 2 (w) = 7(w + w + 1) 7 ~ 2 . 

By C we will denote a constant that depend on the data and thus can be controlled, not necesarily 
arbitrarily small. If the constant depend not only on the data, but also on e, we will denote it by 
C e . Finally, by E we will denote a constant dependent on the data that can be arbitrarily small 
provided that the data is small enough. 
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Since we will usually use the spaces of functions defined on Q, we will omit Q in the notation 
of a space, for example we will denote the space L 2 (Q) by L 2 . The spaces of functions defined 
on the boundary will be denoted by ^(r) etc. 

We do not distinguish between the spaces of vector-valued and scalar-valued functions, for 
example we will write u G W 2 instead of u G (W 2 ) 2 . 



3 A priori estimate for the linearized elliptic system 

In this section we show an e - independent estimate on ||ii £ || W 2 + Hu^H^i, where (u e ,w e ) is a 
solution to (|2.4I) . The first step is an estimate in H 1 x L 2 . Next we eliminate the term divu 
from the second equation applying the Helmholtz decomposition and the properties of the slip 
boundary conditions. Then we derive the higher estimate using interpolation. 



3.1 Estimate in H l x L 2 

In order to prove a priori estimates on H 1 - norm of the velocity and L 2 - norm of the density for 
the system (12.41) let us define a space 

V = {v G ^(QjR 2 ) : vn\ r = 0}. (3.1) 

The estimate is stated in the following lemma. 

Lemma 1. Assume that e, \ \u\\wg and ^\H>\\wi are small enough and f is large enough. Then for 
sufficiently smooth solutions to system < \2.4\) the following estimate is valid 

\\u\\w} + IMUa < C[\\F(u,w)\\ v * + \\G(u,w)\\ L2 + \\B\\ L 2 {r) +E\\w\\ w i]. (3.2) 

where V* is the dual space ofV. 

Before we start the proof, we shall make a remark concerning the term | |iw| that is rather 
unexpected in an energy estimate. Its presence is due to the functions ax(w) and (w + w + 1) 
on the r.h.s. of (12.41) . However, this term does not cause any problems when we apply (13.21) to 
interpolate in the proof of Theorem [21 since it is multiplied by a small constant. 

Proof. The proof is divided into three steps. First we multiply (12.4l h by u and integrate over 
Q. We obtain an estimate on | \u\ |#i in terms of the data and | |u>||l 2 . Then we apply the second 
equation to estimate | \w\ \l 2 and finally combine these estimates to obtain (13.21) . 

Step 1. We multiply (I2.4l) i by u and integrate over Q. Using the boundary conditions (12.41) ^ 
we get 

f Q 2jj~D 2 (u) + i>div 2 udx + J r (f + ^)\u\ 2 da + j Q [ai(w)}S7wudx = 
+ f~ Fudx + f r B(u ■ t) da. 

The boundary term on the l.h.s will be positive provided that / is large enough. Next we integrate 
by parts the last term of the l.h.s of (13.31) . Using (12.41) ? we obtain: 

[ai(w)]Vwudx = — I [a\{w)\ div uwdx — / u w V[ai(w)] dx = 
Jq Jq 
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r ^^)1 Y f \ f 

w 2 n^ da — / w 2 [d Xl a2(w) + (u+uo)Va2{w)]dx — / [a,2(w)]div (u+uq) w 2 dx 

2 2 J Q 2 J Q 

— / [a>2(w)]G(u, w)w dx — e / [a 2 (w)]wAw dx — / uwV[a\{w)]dx. 

JQ JQ JQ 

Since |r out = 1, using 03 .31) and the Korn inequality ( (17.11) . Appendix) we get: 

Cq\\u\\ 2 w i + J ro Ja 2 (w)]w 2 da < 

< / a 2 (w) div (u + Uo)w 2 dx + / [a 2 (w)] Gw dx + / Fudx+ / B(u-r)da 
JQ jq JQ Jr 



"V - 

h h 



+ e / [a 2 (w)]wAw dx + / uw V[ai(u>)] dx + / w 2 [9 xi a 2 (w) + (u + u ) Va 2 (w)] dx . 
jq Jq Jq 

v v ' » „ ' „ ' 

h h h 

(3.4) 

Obviously we have Ii < E\\w\\ L2 . Now we have to deal with the term with Aw. Due to the 
boundary conditions (|2.4lu R we have 

J 3 = e / [a 2 (w)]w)Awdi = -e / [a 2 (w)] | Vu>| 2 dx — e / w V[a 2 (w)] Vw dx. (3.5) 
Jq Jq Jq 

Using Holder inequality we get 

| / wV[a 2 (w)]Vwdx\ < \\V[a 2 (w)]\\ Lp \\wVw\\ Lp * < 
Jq 



< l|V[a 2 («))]|U|VHk \H\ Lq < C\\Vw\\ 2 L2 , 

where q = < +oo and p* = Thus the term with e on the r.h.s of (13.41) will be negative 
provided that | \w\ \ w i will be small enough. Next, 

h < I / uwV[a 1 (w)]dx\ < C||V[ai(«;)]||i, \\u\\ w i IMUa - ^(IMIw/i + I MIL)- 
Jq ~ 2 

The last term of the r.h.s. is the most inconvenient and it must be estimated by - norm of 
w, and this is the reason why this term appears in (13.21) . Fortunately it is multiplied by a small 
constant what will turn out very important in the proof of Theorem [2l We have 



^5 < C r ||a 2 (ty)||wi IMIwj - E W w \\w^ 



Provided that the data is small enough, using the trace theorem to estimate the boundary term 
and the Holder inequality we get 
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n + C L w 2 da < 

" 11 vV$ J I out — /o z:\ 

< C [\\F{u,w)\\ v . + \\G(u,w)\\ L2 + ||S||L a( r)](IHIw5 + IIHUJ + E\\w\\ 2 wl . P " D; 



Step 2. In order to derive (13.21) from (13.61) we need to find a bound on | \w\ \l 2 - From (12.41) 9 
we have 

d xi w = G — (u + u ) ■ Vw — (w + Wo + 1) div u + eAw, 

thus 

w 2 (xi, x 2 ) = w 2 (0, x 2 ) + / 2w w s (s, x 2 ) ds = 



Jo 

/ 2w[G — (w + wq + 1) c?iv u] ds — / 2w 2 (u + uq) ■ Vw 2 ds + 2e / to 2 Aw 2 <ix . 
</o Jo Jo 



Si 5 2 

5i can be estimated directly: 

/ S x < (\\G\\ L2 + C\\u\\ Hl )\\w\\ L2 . (3.7) 

JQ 

It is a little more complicated to estimate S 2 . We have 

pxi px\ px\ 

S 2 = — I [u + u )^d s w 2 (s } x 2 ) ds — / (w + Uo)^d X2 w 2 (s, x 2 ) ds + 2e / wAw. 
Jo Jo Jo 

Now we integrate both components by parts. In the second component we use the fact that the 
integration interval does not depend on x 2 . We get 



S 2 = -(it + u ) w w 2 (x 1 ,x 2 ) + I (u + u )^w 2 (s,x 2 ) ds 

d 







dx 2 Jo 



xi rx\ rx\ 

[u + u )^w 2 (s, x 2 ) ds + / [u + uo) x 2 2 w 2 (s, x 2 ) ds + 2e / wAw 



2:2 

^0 



/•Xl ^ PX1 

(u + Uo)^w 2 (xi, x 2 ) + / if 2 cfiu (w + u )(s, x 2 ) ds — — — / (u + u )^w 2 (s,x 2 ) ds+ 

Jo ox 2 Jo 

XI 

2e / wAw =: Si + S 2 + + St 







2 ~ 2 ~ 2 ' u 2 

The integrals of S\ and S| can be estimated in a direct way: 



/ 1^1, / \S. 
Jo Jo 



2 J , - L • I I. , . I I - 



<E\\w\\l2. (3.8) 



Next, 



Si = I JL [ / u (2 V(s, x 2 ) ds] = / n (2) [ / (u + uo) {2) w 2 (s, x 2 ) ds] da. 



2 Jo 



Now we remind that w = on r in . Moreover, the boundary conditions yields [u + uq)^ 1 = 
on r . Finally, on T out we have = 0. Thus 

S 3 2 = 0. (3.9) 



Finally, 



rl rl pxi rl r 

S^dx = [ / wAw(s, X2) ds dx2\ dx\ = / [ / wAw(x) dx\ dx±, 
Jo Jo Jo Jo Jp xi 

where P Xl := [0, x\] x [0, 1]. We have 

f f f {23J 5 ,6 f 1 

/ wAwdx = — / \Vw\ dx + wVw-nda < I ww xi (xi, 22) dx2, 
Jp xi Jp xi JdP xl Jo 



thus 

-1 rl 



/ S^dx<2e / / ww xi (xi, X2) dx2 dx± = e / d Xl w 2 dx = e / w 2 n^ da. (3.10) 
Jq Jo Jo Jq J r out 

Combining (M . <f3T9l) and (|3A0|) we get 

S 2 = / Si + S 9 2 + Si + S* < E\\w\\ 2 r2 + e I w 2 da. 



2 — / 2 * 2 * 2 ' 2 — N w Nl 2 
q Jq ~ 

Combining this estimate with (13.71) we get: 

HHII2 < C(\\G(u, w)\\l 2 + \\ u \\wi) 2 + E IMU 2 + e / w 2 da, 

J ^ out 

and thus 

\\w\\ 2 L2 <C(\\G(u,w)\\ L2 + \\u\\ w i) 2 + Ce f w 2 da. (3.11) 

Tout 

Step 3. Substituting (13.111) to (UB we get: 

\\u\\ 2 m + I w 2 da < CD(\\u\\ w i + \\w\\ L2 ) + CD 2 + E\\w\\ 2 w i, (3.12) 

^ Tout 

where D = \ \F(u, w)\\v* + \\G(u, w)\\tf+ \ \B\\ L srpy Combining this inequality with (13.111) we 
get 

(\\u\\ w i + \ \w\\ L2 ) 2 + (C - e) / w 2 da < C D (\\u\\ w i + \\w\\ L2 ) + D 2 + E \\w\\ 2 w i, 
thus for e small enough we obtain (13.21) . □ 

3.2 Estimate for 1 1 u \ \ W 2 + 1 1 w \ | ^1 

The following theorem gives an e - independent estimate on ||w e ||v^2 + 1 1 'tt7 e 1 1 where (u e , w e ) 
is a solution to (12.41). 



Theorem 2. Suppose that (u e ,w e ) is a solution to A2.4\) . Then the following estimate is valid 
provided that the data, \ \u\ \ W 2 and \ \w\ \ w i are small enough and f is large enough. 

IMIwj + \\w e \\ w i < C [\\F e (u,w)\\ Lp + \\G t {u,w)\\ w x + \\B\\ w i-i/ V{T) ], (3.13) 

where the constant C depends on the data but does not depend on e. 

The proof will be divided into three lemmas. In the first lemma we eliminate the term div u 
from d2T31)2- 



Lemma 2. Let us define 

H := -(v + 2p) divu € + [ai(w)]w e . (3.14) 
where (u e , w e ) is a solution to A2.4\) and a\{w) is defined in rt2.7D . Then 



\VH\\r < C 



\FJu,w)\\ Lv + ||£|| u/ i-i /p , r . + \\u\ 



*;" 1/p (r) 1 ii"-nw^- 1/p (r)_ 
and w e satisfies the following equation 

[ao(w)]w + w Xl + {u + u ) ■ Vw — eAw = H, (3.16) 

where 

6 = m* + v* + i) + 

v + 2/x 

Proof. Let us rewrite (I2.4l) i as 

d xi u f — fiAu e — (v + fi)Vdivu e + 7Vu> e = F e (u,w) — [ai(w) — 7]Vw e . 

Taking the two dimensional vorticity of (I2.4lh we get 

d Xl a e - fiAa e =rot[F e (u,v) - (a^w) - j)Vw e ) in Q, 

a e = r) + § on T, (3 - 18) 

where a e = rot u e = uf} xi — u^l 2 . The boundary condition (13.18I) ? has been shown in IfTOll in a 
more general case; a simplification of this proof yields (I3.18l) 9. Since our domain is a square, we 
can use the symmetry to deal with corner singularites and apply the standard L p theory of elliptic 
equations ([O) to obtain the estimate 

IKIk p i < ^Ol^(^'^)ll^(Q) + ll( a i(^) _ ^) V ^IUp + ll _ ~K' r ) + ~llw/ p 1 - 1/p (r)]- (3 - 19) 
From the definition of ai(w) (|2.7I) we see that 1 1 (a^w) — 7) | \ Loo can be arbitrarily small provided 



E\\w\\ w i (3.15) 



that ||it7||wi is small enough. Moreover, from the boundary condition (12.4IU we have u e = 
r(u t ■ t) on T, thus (|3.19l) can be rewritten as 

\We\\w^ < C[\\F e (u,w)\\ Lp{Q) + \ \B\\ w i-i /P{v) + \ \u e \\ w i-i, V{v) + E \\w e \\ w i\. (3.20) 
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Now we apply the Helmholtz decomposition in of u e (see Appendix, (17.21) ): 

« £ = V0 + V- L A (3.21) 

For simplicity we omit the index e in the notation of <p and A. We have n ■ V L A = r ■ X7A = -j^A, 
thus the condition n ■ V ± A\y = yields A\ r = const. Moreover, 

rot u = rot{V<p + V^A) = rot V^A = AA. 

We see that A is a solution to the following boundary value problem: 

f AA = a e E W}(Q), 
\ A\r = const. 

Applying again the elliptic theory we get 

||-4||w3(Q) < ll«llwj(Q) - C {\\ f ^ u ^)\\l p (Q) + \ \ B \\w^- 1/p (r) + IKHw^cr)} ■ ( ^ 2T > 
Substituting the Helmholtz decomposition (17.21) to (|2.4l) i we get 

d Xl {V(j) + \7 ± A) -/iA(V0 + V^) - (v + /j,)Vdiv(V(f) + V ± A) + [ai(w)]Vw = F e {u,w), 
but divVcj) = A(f) and thus 

-{y + 2//) VA0 + V(K (w)]w) = 
= F(u, w) + fj l AV ± A + (/i + u)VdivV ± A - d Xl V ± A + d Xl V(f) + u;V[ai(tD)] =: F, 

(3.23) 

what can be rewritten as: 

V {-{u + 2/i)A0 + yai(w)w) = F. 

We have A0 = div u, thus F = VH where H is defined in (I3.14I) . From (I3.23I ) we have 

\\F\\ Lp < C[\\F e (u,w)\\ LP + \\A\\ w s + ||V 2 0||lJ + ||V[a (^)]|UJ|«; e ||oo < 
< C[\\F e (u,w)\\ LP + \\A\\ W 3 + M\ W 2]+E\\w e \\ w i 

and from (I3.22I) and (I7.3I) we get (I3.15I) . The proof is thus completed. □ 

In the next lemma we will use the equation (I3.16I ) to estimate | \w\ \ w i in the terms of functions 

H and G(u, w). 

Lemma 3. Under the assumptions of Theorem^the following estimate is valid: 

<C[\\H\\ w , + \\H\ T J\ Lp[Tm) l (3.24) 

where 

H = ^— + G(u,w). (3.25) 
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Proof. Throughout the proof we will omit the index e denoting w e by w. The proof will be 
divided into four steps. First we estimate | \w\\i , then \\w Xi \\l and||w X2 ||z, and finally combine 
these estimates. 

Step 1. Multiplying (13.161) by \w\ p ~ 2 w and integrating over Q we get: 

ao(w) \w\ p + / \w\ p ~ 2 w w Xl + / (u + u ) ■ Vw\w\ p ~ 2 w — e / Aw\w\ p ~ 2 w = / H\w\ p ~ 2 w . 
Jq Jq Jq Jq 



il q if if if 

(3.26) 

We have 

if = - [ (u + u ) ■ V\w\ p dx = — / \w\ p div (u + u ) dx H — / u^\w\ p da. 
PJq PJq PJr out 



Next, 



if = 1 f d Xl \w\ P = - I \w\ p n^ = - f 
PJq PJr P Jr 

Combining the last two equations we get 



\w 



v 



-(I 2 + l!)<E\\w\\l p -C f (l + u^)\w\ p da. 



The boundary term is positive due to the assumption of smallness of u . The term with Aw: 

dw 



I{ = e [ Vw-V(\w\ p ~ 2 w)-e [ \w\ p ' 2 
Jq Jr 



dn 

The boundary term vanishes due to the conditions (12.41k fi and the first term of the r.h.s is equal 
to 



C f{\Vw\ p - 2 w 2 + \w\ p ~ 2 \Vw\ 2 )dx > 0. 
Jq 

The r.h.s of (13.261) can be estimated directly: 

i!= [ H\w\ p - 2 wdx<\\H\\ Lp [ (M (p ~ 1)p y /P *HI#IUMIC- 

Jq Jq 

The smallness of w and wo in implies that ao(w) > C > 0, thus combining the above 
estimates we get C \\w\\ p Lp < II-^IUpIMIl" 1 + E \ \w\\ p L , thus 

\\w\\ Lp <C\\H\\ Lp . (3.27) 

Step 2. In order to estimate w X2 we differentiate (|3.16l) with respect to x%, multiply it by 

\w X2 \ p ~ 2 w X2 and integrate over Q. We get 



Q 



a (w)]\w X2 \ p + / [a (w)} X2 w\w X2 \ p w X2 + / \w X2 \ p 2 w X2 w X2Xl 



^2 ^2 ^2 
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/ {{u + u ) X2 -Vw)\w X2 \ p 2 w X2 + / {{u + u )-V 
Jq Jq 



Aw X2 \w X2 \ p 2 w X2 = / H X2 \w X2 \ p 2 w X2 . 



4 4 



We have 



l\ = - / <9 xl |u> X2 | p <ir = — / |u; X3 | p g?ct H — / \w X2 \ p da, 



V Jq ' V Jr in ' V Jr a 

but the condition w = on Y in implies w X2 = on T in , thus 



P Jv out 

Obviously we have If < E \ \ Vw\ \ P L . Next, 



Il=- ( M p da. (3.28) 



(!)„(!) 



71 v 'ItWa-j, 



if = — / div(u + u ) \w X2 \ p H — / Mq 
P P Jr out 

Combining this equation with (13.281) we get 

II + l\ = -1/p I div{u + u Q )\w X2 \ p + - / (1 + u^KJ'tZo-, 
■/ Q P /rw 

The boundary term is nonnegative due to the smallness of u . 
The last part of the l.h.s: 

il = -e J Aw X2 \w X2 \ p w X2 = e J Vw X2 ■ V(\w X2 \ p ~ 2 w X2 ) + e J ^^-\w X2 \ p ~ 2 w X2 da. 

The first term equals J Q (p — 1) \w X2 \ p ~ 2 \ Vw X2 \ 2 dx > and the boundary term vanishes due to 
the boundary condition (|2.4b /i ^. Using the definition of a (w) (12.71) we get 

/ [a {w)) Xi w\w Xi \ p ~ 2 w Xi <C\\(iD + Wo)^\\L p \\wx i \\j^ 1 \\w\\w^ <£|Ml^ p i, (3.29) 

J Q 

thus 1\ < E Hwll^!. In order to estimate the r.h.s we use the definition of H and the Holder 
inequality. We get 



l\ = I / H x .\w Xi \ p 2 w Xi dx\ < C \\H\\ w i\\w 2 
Jq 



m-l 



(3.30) 



The important fact that we could write H instead of H on the r.h.s easily results from the defini- 
tion of H (|3.17l) . Combining the above estimates we get 

ll^lli^C^HV^II^ + Cllifll^ll^ll^ 1 ]. (3.31) 
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Step 3. In order to estimate w Xl we differentiate (12.41) with respect to x\ and multiply by 



a (w)\w Xl \ p + / [a (w)] xl \w Xl \ p 2 w Xl + / w XlXl \w Xl \ p 2 w Xl + / (u + u ) ■ Vw Xl \w Xl \ p 1 

Q JQ S 



4 4 
- / eAw Xl \w Xl \ p ~ 2 w Xl = / H Xl \w Xl \ p ~ 2 w Xl - (u + u ) Xl ■ Vw\w Xl \ p ~ 2 w Xl . 

JQ JO 



[5 t6 
J 3 J 3 



We have 

ll = - I d T Aw r A p dx = — / Im™ | p dcr. 



l 3 — ~ / "X1\<" X1 — I UA^ 



Next, 

/* /* d iv 

-if = e y Vu^ • Vd^JP" 2 ^.,) dx-ej— n 

The first term is nonnegative and the boundary term reduces to: 



— \w xi \ p 2 w Xl da. 



e / w XlX ,\w x ,\ p 2 w Xl da. (3.32) 



r, 



Note that on Y in equation (13.161) takes the form: 

(1 + u 1 + uQw Xl - ew XlXl = H\r in - 
Thus (13.321) can be rewritten as 

/ [(1 + u 1 + ul)\w Xl \ p - H\w Xl \ p - 2 w Xl } da. 

J Tin 

Finally, 

In = — / diviu + uo)\w xl \ p dx / u\\w xi \ p da. 

PJq P Jr in 

Combining the above results we get 

C [ \w Xl \ p dx+ [ (1 - u 1 - -)\w xi \ p da < 
Jq Jr in P 

<- I div(u + u )\w xi \ p dx + I H Xl \w Xl \ p ~ 2 w xl dx - f [a (w)} Xl w\w Xl \ p ~ 2 w Xl dx 
P Jq Jq Jq 

- / (u + u Q ) xl ■ Vw\w xi \ p ' 2 w xi dx + / H\w xi \ p ~ 2 w xi da, 
Jq Jr in 
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thus using (13.301) and (13.291) we obtain 

(c - f)l Kl l^ W ) + (1-^^) IK, 1 1^)< 
< C ll^llwpi I l^x I li p( Q) + #1 1 VH li p(Q) + II^IU^) I \w Xl 1 1^^) + e I h 

Step 4. Combining (13331) and (I33TT) we get 



(3.33) 



|P +11?/) Il p < 



< C [(HffHwj + IMU IIV^H^ 1 + ll^llMr^lKJI^Vj. 
Combining this estimate with (|3.27l) we get 

\H\l v + \\Vw\\l p + \\w Xl \\ Lp( r in ) < C(\\H\\ Lp + \\H\\ w i + ||tf|| M r in )). (3.34) 

Due to (13.301) we have | \H\ \ w 1 instead of | \H\ \ w 1 on the r.h.s. and the proof of (13.241) is almost 
complete. Now it is enough to note that due to the smallness of w and wq in Wp we have 

\\H\\ Lp < C\\H\\ Lp and \\H\ Tin \\ Lpi r\ in) < C\\H\ Tin \\ Lp{T \ in) , 

thus (I3J41) is indeed <f3T24|) □. 

In order to complete the proof of Theorem[2] we have to estimate H. We will make use of the 
interpolation inequalities (Lemma [TT] in the Appendix). 

Lemma 4. Under the assumptions of Theorem\2\ V<5 > the following estimate is valid 

\\H\\w$(Q) + H-frlrJUpOV) < 

< 5\\u\\ W 2 + C(8)[\\F e (u,w)\\ Lp + \\G e (u,w)\\ wi + \\B\\ w i-i /P(n + E\\w\\ W i], ^ } 

where H is defined in ( 13.251) . 

Proof. For simplicity let us denote F := F e (u, w) and G := G e (u, w). Applying the interpo- 
lation inequality (17.41) to the term | \u\ \wj in (13.151) we get: 

Wp" 1/p {v) + H M Hw^ 1/p (r) + ^il 1^1 lw-2 + C(^i)||w||hi] 
+ E \\ w \\wj- 

In order to estimate we need to apply the interpolation inequality (17.41) and then the 

energy estimate (13.21) . We get 

\\H\\ Lp < 5 2 \\VH\\ Lp + C(S 2 )(\\F\\ L2 + ||G|| La + ||£|| i2(r) + E\\w e \\ W i). 

Combining the above estimates we get 

||V#|| Lp + \\H\\ Lp < 
< (l + 8 2 )\\VH\\ Lp + C(5 2 ) l\\F\\ P L * + \\G\\v + \\B\\ L2{r) + E\\w £ \\ w i] < 

< S 3 \\u e \\ W 2 + C(S 3 )[\\F\\ Lp + \ \G\\ w i + \ \B\\ w i- 1/P{r) + llwell^i-i/p^j +E\\w e \\ w i]. 

(3.36) 
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Using the trace theorem, (17.41) and (13.21) we estimate the boundary term \ \u\ 



(iy 



Ur 



i w ^/p [T) < $4\ue\\w* + C(5 4 )[\\F\\ L2 + \\G\\ L2 + \\B\\ L 2 {V) + E\\w £ \\ w i\. (3.37) 



In order to complete the proof it is enough to estimate H\ r . n . By the trace theorem we have 
\\divu\\ Lp{rtn) < C{r)\\divu\\ w i/ P +r {Tm) Vr > 0. Thus, since w\ Tm = 0, applying ([73]) and 
(Q we get 

||^|ri„IUj,(r in ) < \ \divu e \ Ttn \\ Lp (r m ) + | \G\ v . n \ \L p (v m ) < 
< C r \\divu t \\ w i/ P +r {Q) + \\G\\ w i( Q ) < S 5 \\u e \\ W 2 + C(8 5 )\\u e \\ Lp + \\G\\ w i < 

< 5 5 \\u € \\ W 2 + C(6 5 )[5 6 \\u\\ W 2 + C(S e )\\u\\ H i) + \\G\\ w i < 
< 8 7 \\u e \\ W 2 + C(8 7 )[\\F\\ L2 + \\G\\ L2 + \\B\\ L 2 (r) +E\\w e \\ w i]. 



(3.38) 



Since we control the smallness of <5 5 , we also control C(5 5 ) and thus we can choose 5 6 to make 
87 = 5 5 C(5 6 ) as small as we want. Next, substituting (13.371) to (13.361) with (??) we get (13.351) 
with S arbitrarily small since 5\ . . . 67 can be arbitrarily small □. 

We are now ready to complete 

Proof of Theorem\2\ Let us fix rj > 0. Provided that w and w are small enough, combining 
(1X241) and (13351) we get 

IKIki(Q) + IK.ailUpCr*,) < n oqn 

< nWuWwg + C(rj)[\\F e (u,w)\\ Lp + \\G e (u,w)\\w$ + \\B\\ w i-i /P{r) ]. ^> 



The theory of elliptic equations applied to (I2.4l) i yields 

WueWw* < \\F e (u,w)\\ LP + \\w t \\ w i. (3.40) 
Combining this estimate with (I3.39I ) we get 

\\ u e\\w* + IKHiyi < T]\\u e \\ W 2 + C v [\\F £ (u,w)\\ Lp + \\G e (u,w)\\ w i + \\B\\ w ^i /P{r) ]. 
Choosing for example rj = | we get (13 .13b . □ 

4 Solution of the linear system 

In this section we will show that the operator T e is well defined. We have to show that the system 
(12.41) has a unique solution (u, w) G W? x for (u, w) G x small enough. The 
necessary result is stated in the following 

Theorem 3. Assume that \ \u\\ W 2 + | is small enough. Then the system ( 12.41) has a unique 
solution (u e , w e ) G Wp x W£ and the estimate < \3.13\) holds. 
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We shall make here one remark concerning the above theorem. The fact that (u e ,w e 



W„ x Wp is a consequence of the ellipticity of the system (12.41) . but the estimate on ||u>||w"2 



depends on e. What will be crucial for us is that (13.131) does not depend on e. 

The system (12.41) has variable coefficients thus its solution is not straightforward. In order 
to proove Theorem |3] we will apply the Leray-Schauder fixed-point theorem. Given (u, w) £ 
W 2 x Wl we define an operator S e (a ^ : W 2 x W 2 -> W 2 x W%\ 
(u, w) = St- ^ (u, w) -<=>- (u, w) is a solution to 

d Xl u - fiAu — {v + n)Vdivu + 7Vw = FL^Ju, w) 
div u + d Xl w - eAw = G^-^ (u, w) 
n ■ 2fiD(u) ■ r + / u - t = B 
n ■ u = 
w = 

dw Q 

dn 

where 



in 


Q. 


in 


Q, 


on 


r. 


on 


r, 


on 


r. 

L mi 


on 


r \ r 

x \ 1 mi 



(4.1) 



F (u,w)( u ' w ) = -(ai(«0 -7)Vw + F e (M,w), 
= — (w + wo) rf^M - (w + u Q ) ■ Vw + G e (u,w). (4.2) 

We have to show that Sf~ ^ is well defined and verify that it satisfies the assumptions of the 
Leray-Schauder theorem. The reason to consider Sf^ ^ on W 2 x W£ instead of W 2 x is that 
it is straightforward to show its complete continuity. 

4.1 Solution of the system with constant coefficients 

In this section we show that the operator SL is well defined. Thus we have to show that the 
system 



d Xl u — fiAu — 


{y + ji)Vdivu + •y'Vw = F 


in 


Q. 


divu + d xi w - 


- eAw = G 


in 


Q. 


n ■ 2/iD(u) ■ t 


+ f u ■ T = B 


on 


r, 


n ■ u = 




on 


r, 


w = 




on 


r 


dw Q 

dn 




on 


r \ r 



(4.3) 



where F,G £ W}(Q) are given functions, has a unique solution (u, w) £ W 7 ^ x Wi 2 . We start 
with showing existence of a weak solution to the system (14.31) . Let us recall the definition of 
space V (13.11) and introduce another functional space W = {w £ H l (Q) : w\r in = 0}. Consider 
a bilinear form on (V x W) 2 : 
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—7 / wdivvdx + 'y / r/divudx + 'y / rjd xi w dx 4 / Vw-Vr]dx 
Jq Jq Jq Jq 

and a linear form on (V x W): 

F(v,7)) = / F-vdx+ / B(v-r)dx+ / Grjdx. 
Jq Jr Jq 

By a weak solution to the system (14.31) we mean a couple (u,w) E V x W satisfying 

B[(u,w),(v,r])} = F(v,r]) V(v, 77) e V x W. (4.4) 
Using the definition of V and we can easily verify that 

H[(u,w),(u,w)]> / 2fiD 2 (u) 4 vdiv 2 udx 4 e / | Vu>| 2 c?x > C e [\ \u\ Ih^q) + I \ w \ Ih 1 ^)], 



thus existence of the weak solution to (14.31) easily follows from the Lax-Milgram lemma. Using 
standard techniques we show that the weak solution belongs to W 2 (Q) x W 2 (Q) and 

IMIw? + \\w\\w2 < Ce[\\F\\ Lp 4 \\G\\ w i 4 \\B\\ w i-i/ P{r) ]. 

4.2 Complete continuity of Sf- - A 

In this section we show that SL ^ is continuous and compact. Since it is a linear operator, it is 
enough to show its compactness, and this is quite obvious due to elliptic regularity of the system 
(14.31) . Namely, if we take a sequence (u n , w n ) bounded in W 2 x W 2 , then the sequence 

(F(u,w) {u n , w n ) , Gfaw) {u n , w n )) 

is bounded in x W^. Thus the sequence (u n , w n ) = SL iD \(u n , w n ) is bounded in W£ x W£ 
(the bound on ||w|| W 3 depends on e, but at this stage e is fixed, so it does not matter). The compact 
imbedding theorem implies that (u n , w n ) has a subsequence that converges in W 2 (Q) x W% (Q). 
Thus SL ^ is compact. 

4.3 Leray-Schauder a priori bounds 

Next we have to show a A - independent a priori estimate on solutions to the equations (u\, w\) = 

^ S (u,w)( U ^ W *)> that read 



(4.5) 



for A G [0, 1]. Actually we should write (u\, w\), but we will omit e as it should not lead to any 
misunderstanding. The result is stated in the following 



d Xl u x - nAu x - 


- (y 4 i^Vdivux 4 


- 7 Vw A = \F{ a{d) (u x ,w x ) 


in 


Q. 


div u\ 4 d Xl W\ 


— eAw\ = XG^y 


[ux,w x ) 


in 


Q, 


n ■ 2{iD(u\) ■ t 


+ fu x -r = B 




on 


r, 


n ■ u\ = 






on 


r, 


w x = 






on 


r 


dw x _ q 

dn 






on 


r\ r 

1 \ x in) 
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Lemma 5. Let (u\,w\) = \S e { ~ ^(u\,w\), then 

IMIw* + ||wa|Iw* < c e [\\F(u,w)\\l p + \ \G(u,w)\\ w i + \ \B\\ w i-i /P{r) }. (4.6) 

Proof. The proof is very similar to the proof theorem [2l First we repeat the proof of Lemma 
Q] obtaining the A-independent energy estimate 

\\u\\\m + IKIk < C0I F (^)( m a,w a )|| L2 + \\Gl u>iD) (ux,w x )\\ L2 + ||-B|| i2(r) ] +E\\w\\ w i. 

(4.7) 

Next we take the vorticity of (14.51) : 

d Xl a x - |uAa A = rot (AF^^.wa)) in Q, 
a\ = -fax-T) + f on T, 

where a x = rot u x . Thus 

HatAllwf < c {\\ F (u,v)( u \i w x)\\l p (Q) + \\ B \\ W ^- 1/P (r) + I \ Ux \ \wp~ 1/p (r) } ' 

Now let u x = V0a + A A . Substituting this decomposition to (14.51) we get 

-(v + 2/i)VA</> A + V( 7 ^a) = A Ffc^ux, w x ) + - d Xl A{ + <9 Xl V0 A - =: F A , 

what can be rewritten as: V (— [y + 2/i)A0 A + 7W A ) = F A . We denote as previously 

{—{v + 2fji)dwux + [ai(w))wx) = H x . 

Combining this identity with (14.51) 9 we get an analog of (13.161) : 

(x(w)wx + w x , xi + X(u + u ) ■ Vw x - eAwx = Hx, (4.8) 

where (x(w) = [1 + X(w + Wq)\ and Hx = 1+ ^ 2 ^ u ' o) i/ A + AG. Now we can repeat step by 
step the proof of Theorem [2] obtaining the estimate 

IMIwj(Q) < ^II^aIIw^ + Cr?[II^^,-a>) (-"a, ^a)| U y -h I iG^c^) (-"a, -wa)Hwi "I- ll^ilw^-^cr)] ( 4 - 9) 

for each rj > 0. The estimates for | \u\\ \wjj an d I \ w x\ I w| now easily result from the system (I4.5I ). 
Namely, applying the standard elliptic theory to (14.5lh we obtain an estimate 



|«a||w/ 2 < c 



(4.10) 



M\w$ + \\F(u,w)( u ^ w \)\\l p 

that does not depend on A. Next, from (|4.5I> 9 we get an elliptic estimate 

IMIw/2 < C e (\\wx\\wi + \\u\\\wg + \\ g Iu,u,)( u \i w x)\\l p )- (4.11) 
Combining (l4~9l) . (14.101) and (14.111) we get 



[H^H^-^cr) + W F (u,w)( u ^ w x)\\l p + \\G c (iitiB) (ux,wx)\\ Lp ], (4.12) 
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but from the definition of Ft ^ and G e ,- ^ we have 

II F (W m a,wa)||l p + \ \G\u,w)( u ^ w x)\\l p < 

E (\\u\\\w* + ||wa||w«) + \ \Fe(u,w)\\ Lp + \\G £ (u,w)\\ Lp 

and thus (I4T21) yields (1431) . □ 

Now we are ready to complete 
Proof of theorem |3l We have shown that the operator SL ^ satisfies the assumptions of the 
Leray-Schauder theorem. Thus there exists a fixed point (u e , w e ) = S^-^(u e ,w e ). The fixed 
point is a solution to (12.41) . Its uniqueness follows directly from the estimate (13.131) . □ 

We have shown the existence of a unique solution to the system (|2.4I) under some smallness 
assumptions on u and w. Thus we define the domain V of the operator T: 

V = {(u,w) G W*(Q) x Wp(Q) : Theorem[3]holds for (u, w)}. (4.13) 



5 Solution of the regularized system 

In this section we show existence of a solution to an e-elliptic regularization to the system (12.31) . 
The result is stated in the following 

Theorem 4. Assume that the data and e > are small enough and f is large enough. Then there 
exists a fixed point (u*, w*) = T t (u*, w*) and 

\K\\w* + \M\\w}<M, (5.1) 

where M depends on the data but does not depend on e and can be arbitraily small provided that 
the data is small enough. 

In order to prove the theorem we apply the Schauder fixed point theorem to the operator T € 
defined in ((231) . 

Lemma 6. Assume that uq and wq are small enough. Then T e (B) C B for some ball B C 
W 2 p {Q) x W}{Q). 

Proof From the definition of F e (u, w) and G e (u, w) we have 

\\F e (u,w)\\ w i + \ \G e (u,w)\\ w i <E + (\\u\\ W 2 + H^H^i) 2 . (5.2) 
Thus we can rewrite the estimate (13.131) as 

IH|w| + IKIIw* < C [D + (\\u\\ W 2 + ||w||wi) 2 ], (5.3) 

where D can be arbitrarily small provided that ||« Q | \ W 2 and | |w | \w^ are small enough. In (13.131) 
we only need an estimate on ||F e ('u,'iD)||L that hold also for F{u,w), but we will need the 
estimate in to show the compactness of T e and this is the reason why we introduce the 
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regularization F e . Let us assume that the data is small enough to ensure D < -^p, where C and 
D are the constants from (15.31) . Assume further that | \u\ \ W 2 + \\w\ \ w i < yD. Then from (15.31) 
we get 

I \ u e\ \wg + 1 1 We | IWp 1 < 2C D < \f~b. 

Thus T e (B) C B where B = B(0, VD) C W^ 2 (Q) x W^Q) □. 

In the next lemma we show that T e is a continuous operator on V, where V is defined in 
(14.131) . The proof applies the estimate (13.131) which requires some smallness assumption, but this 
assumption is also included in the definition of V and therefore we can prove the continuity on 
the whole V. 

Lemma 7. T e is a continuous operator on T>. 

Proof Let us have (ui, Wi) = T(ux,Wi) and (u 2 , w 2 ) = T(u 2 , w 2 ), then the functions u 1 —u 2 
and w\ — w 2 satisfies the equations 

d Xl («i - u 2 ) - fiA(ui - u 2 ) — (v + jj,)Vdiv (m - u 2 ) + 7(^1 + w + 1) 7_1 V(w l - w 2 ) = 
= F e (ui, Wi) - F € (u 2 , w 2 ) - 7[(tDi + w + l) 7-1 - (w 2 + w + l) 7 " 1 ] Vw 2 

and 

(wi + w + 1) div (u\ - u 2 ) + d Xl (w\ - w 2 ) + (u\ + u ) ■ V(wi - w 2 ) - eA(wi - w 2 ) = 

= G e (ui,wx) - G e (u 2 ,w 2 ) - (wi - w 2 ) divu 2 - (u\ - u 2 )) ■ Vu> 2 , 

supplied with boundary conditions 

n ■ 2jiD{ui — u 2 ) ■ t + / (ui — u 2 ) ■ t = on T, 
n ■ (ui - u 2 ) = on T, 

w\ - w 2 = on T in , 

d{w\—W2) 



dn 



on r\r in . 



If (ui,Wi), (u 2 ,w 2 ) E V then the system on {u\ — u 2 ,W\ — w 2 ) satisfies the assumptions of 
Theorem [2] and thus (13.131 ) yields 



||«1 - U 2 \\ W 2 + \ \W! - W 2 \\ W i < 

||F e (Mi,Wi) - F e (u 2 ,w 2 )\\ Lp + WGeiuuW!) - G e (u 2 , w 2 ) I \ W l 
+\\(wi + w + 1) 7_1 - (wi + w + l) 7-1 Vw 2 |U P + \\(w! - w 2 ) divu 2 \\ w i 

+||(«i - u 2 ) ■ Vw 2 \\ w i. 

From the definition of F e (u, w) and G e (u, w) we directly get 

||F e (ui,tBi) - i^ e (it2,«>2)|lx» + ||G e (TZi,u;i) - G e (u 2 ,w 2 )\\ w i 

+\\(wi + wo + 1) 7_1 - (wi + wo + 1) 7_1 Vw 2 ||l p + \\(u>i - w 2 ) divu 2 \\ w i < 
< C (IJuiHwi, HwiH^i, [jMallwj, [|«>2||wj) - U 2 \\ W 2 + \\wi - W 2 \\ w i\. 
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(5.5) 



In order to estimate the last term of the r.h.s. of (15.51) we have to use higher norm of W2: 

\\(ux - u 2 ) ■ VwiWwi < C{\\w 2 \\w$) IK - u 2 \\ w $- 
Since on this level e is fixed, we can use the elliptic regularity of the system (12.41) that yields 

||w2||h/2 < C e \\F e (u 2 ,w 2 )\\ Lp + \\G e (u 2 ,w 2 )\\ w i + ||-B|U P (r)- 
Combining the above estimates we get from (15.51) : 

IK - u 2 \\ w i + ||wi - w 2 \\ w i < Ce [IK - u 2 \\ W 2 + \\w 1 - w 2 \\ w i\, (5.6) 

what completes the proof of continuity of T t □. 

Now we need to proove that T e is a compact operator. The key is in the following lemma 

Lemma 8. Let us have (u, w) = T £ (u, w). Then (u, w) e Wp(Q) x Wp(Q) and 

\H\w3 + lklk| < C e [\\u\\ W 2 + \ \w\\ w i +E]. (5.7) 

Proof. If (u, w) = T € (u, w) then in particular w satisfies 

—Aw = G e (u, w) — d Xl w — (u + uo) ■ Vw — (w + wo + 1) div u. 

Thus by (|3.13l) we have 

\\w\\wi < C e [I \G e (u,w) I \ Lp + C{\\u\\ W 2 + llwllwx)] < 
< C e [\\F e (u, w)\\l p + \\G € (u, w)\\ w i + \ l-BH^i-i/,^] • (5.8) 

Next, u satisfies the equation 

—//Aw — (is + n)Vdivu = F e (u, w) — d Xl u — ^(w + wo + l) 7 ~ 1 Vw , 

what yields 

IMIwj < C [\\ F e(u,w)\\ W l + \ \wWwg] < 

EH 

< C e [\\F e (u, w)\\ w x + \\G e (u,w)\\ w i + \\B\\ w i-i /P(p) ]. (5.9) 

Now, from (I5T21) we get (15771) . □ 

With Lemma[8]the compactness of T e is a straightforward consequence of the compact imbed- 
ding theorem. Namely, if we take a sequence (u n ,w n ) that is bounded in Wp(Q) x Wp(Q) 
and consider (u n ,w n ) = T € (u n ,w n ), then from (|5.7I) the sequence (u n ,w n ) is bounded in 
Wp(Q) x Wp(Q). Thus the compact imbedding theorem implies the existence of a subsequence 
(u nk ,w nk ) that converges in Wp (Q) x Wp(Q), what means that T e is compact. 

Proof of theorem 31 The theorem results directly from the Schauder fixed point theorem for 
the operator T e . □ 
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6 Proof of Theorem U 

In this section we prove our main result, Theorem [Q passing to the limit with e in (12.41) . The 
proof will be divided into two steps: the proof of existence of the solution and the proof of its 
uniqueness. These steps are quite separated since in order to prove uniqueness will will go back 
to the original system (|l.ll) and modify the proof of the estimate (13.21) . 

Step 1: Existence. Consider a decreasing sequence e n — ► 0. If e\ is small enough that 
Theorem[4]holds (what we can assume without loss of generality), then for each n E N Theorem 
|4] gives a solution (u £n , w 6n ) to an e n - elliptic regularization to (|2.3I) . 

By (15.11) the sequence (u 6n , w en ) is uniformly bounded in Wp x Wp. The compact imbedding 
theorem implies that there exists a couple (u, w) E x W^ such that (up to a subsequence) 

wl Wp 
u tn u and w en w. (6.1) 

From the definition of F e and G e we easily get 

F e (u e ,w e ) F(u,w) and G € (u e ,w e ) G(u,w). (6.2) 
We have to show that (u, w) satisfies the system (12.31) . Clearly we have 



Au en Am, Vdivu En Vdivu, 



(6.3) 



Thus it remains to show convergence in nonlinear terms, but this is also straightforward. We 
have V</> E L q : 

/ </>(u> e + Wo + l) 7 ~ 1 Vw e dr = 
Jq 

I 4>[(w e + w + ly- 1 -(w + w + ly- 1 ] vw e + 1 4>(w + w + \y- l vw e dx 

Jq Jq 

Since <j>(w + w + 1) 7_1 G L g , the second term converges to j Q 4>(w + w + l) 7_1 Vw dx. The 
first term 

/ (f)[(w € + w + I) 1 ' 1 - (w + wo + iy- 1 ] \7w e dx < 

Jq 

< Miwe + w + iy~ l ~(w + w + iy- l ]\\ Lq \\w t \\ W i 0, 

since by the compact imbedding theorem w e -4 w V 1 < g < +oo. Thus 

/ 4>(w e + w + I)"*' 1 Vw e dx -> / 4>(w + w + ly^Vwdx. (6.4) 
>/ q Jq 

Similarily we can show that 

(w e + wq + 1) dit> it e + (-u e + -uq) • Vty e (w + wq + 1) div u + (u + uq) ■ Vw. (6.5) 
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From (16.21) . (16.31) . (16.41) and (16.51) we see that (u, w) satisfies (12.3lh p a.e. in Q. The trace theorem 
implies that 

i L p (Ti n ) . . L p (r) . , . L p (r in ) sr\ 

We\ Jin w\ Tin , u\ r -» u| r , D(u) u r . 

Thus it satisfies (|2.3b ? /i a.e. on T and w satisfies (12.31k a.e. on T in . Now take t> = u + u + v 
and p = io + w + p, where u and w are extensions to the boundary data defined in (12.11) and 
(v, p) = ([1, 0], 1) is the constant solution. Then (v, p) satisfies the system (11.11) . 
In order to show the estimate (11.21) we repeat the proof of Theorem [2] obtaining 

\H\w$ + IHIwj ^ C [\\ f ( u i w )\\l p + \ \G(u,w)\\ w i + \ \B\\ w i-i /P{r) ]. (6.7) 

We have 

\\F(u,w)\\ Lp + \\G(u,w)\\ w i <D+ (\\u\\ W 2 + \\w\\ w i) 2 , (6.8) 

where D can be arbitrarily small provided that the data is small enough. From (|6.7I) and (16.81) we 
conclude (|1.2I) . 

Step 2: Uniqueness. In order to prove the uniqueness of the solution in a class of small 
perturbations to of the constant flow (v, p) consider (fi, pi) and (t> 2 , p-i) both being solutions to 
(11.11) satisfying the estimate (|1.2I) . We will apply the ideas of the proof of the energy estimate 
(13.21) in order to show that 

IK -^Hhi + \\p\ -P2||| 2 = 0. (6.9) 

For simplicity let us denote the differences u := v± — v% and w := pi — pi- We will follow 
the notation of constants introduced before, namely E shall denote a constant dependent on the 
data that can be arbitrarily small provided that the data is small enough, whereas C will denote 
a constant dependent on the data that is controlled, but not necessarily small. In order to show 
(16.91) it is enough to prove that 

\\u\\ H i < E\\w\\ L2 (6.10) 

and 

I \w\ \l 2 < C\ \u\ \ H i. (6.1 1) 

If we substract the equations on (vi, pi) and (d 2 , p 2 ) there appears a term pj — p],. We will use 
the fact that pi, p 2 ~ 1 =^ p\ — p\ ~ 7(pi — p 2 ), more precisely, we can write: 

P\-Pl = (Pi - P%) t l[tpi + (1 - t)p 2 ] 7 - 1 ctt 



o 



and we have L ( ~ 7. Now we easily verify that the difference (u, w) satisfies the system 

w v 2 ■ Vt> 2 + pi u • Vf 2 + Pi f 1 • Vw — pA-u — (p + v)Vdiv u + J 7 Vw = 0, 
pi c?ii> m + w div f 2 + m • Vp 2 + f 1 • Vw = 0, 

n-2pD(M)-r| r = 0, (6.12) 
n ■ u\y = 0, 

w\r. = 0. 



24 



We modify the proof of (13.21) . multiplying (I6.12l) i by p\ u and integrate over Q (the reason why 
we take p\ u instead of u will be explained soon). We get 

/ (2pD 2 (w) + u Pl div 2 u)dx+ / [(p x - l)D(u) : Vu + D (u) : (u £g> Vpi)] dx 
Jq Jo 



— J 7 / wp 1 divudx+ / p\ f u da 
Jq Jy 

I wuVpi+ / p\u 2 -Vv 2 dx + / w u> pi f 2 • Vf 2 cte + / p^ (v i • Vu) • udx = 0. 
Jq Jq Jq Jq 




We have \Ii\ + \I 2 \ + + I/4I < E (||w||^-i + ||Hli 2 )- Now let us split J 5 into two parts: 

2J 5 = I {plv? -l)d Xl \u? + pW? d X2 \u\ 2 dx+ [ d Xl \u\ 2 dx. 
Jq Jq 



We have < E\\u\\ 2 H1 and I 2 = L lu^n^ da. The last term can be integrated by parts and 
combined with the boundary term involving friction. Thus applying the Korn inequality (17.11) we 
get 

^IMIiJi+ (pi f + n )\ u \ 2 da — 1^ w divudx < E\\u\\ 2 H i. 
Jv Jq 

For the friction coefficient / large enough the boundary term will be positive and thus 

IMI/p — C I wpidivudx. (6.13) 
Jq 

The reason why we multiplied (|6.12lh by p\ u is that now we have this function on the r.h.s of 
(16.91 ) instead of divu. In order to derive (16.101) from (16.131) we express pi div u in terms of w 
using the equation (16.121b . Thus we can rewrite (16.131) as 

IM|#i<— / w 2 divv 2 dx — / wvi-Vwdx— / wu-Vp 2 dx. (6.14) 
Jq Jq Jq 



h h h 



Obviously |/ 6 | < and, since p > 2, we have \I S \ < ||Vp2||z, p IMU2 ll M IU<, f° r some 

q < 00. Thus from the imbedding theorem we get |/ 8 | < -E'(||tu||i 2 + IM|#i). Integrating by 
parts in I 7 and using the boundary conditions we get 

—2/ 5 = w 2 div Vi dx — / v 1 da. 

JQ JY ou t 



25 



The boundary term is positive since v± ~ 1, thus — J 5 < C \ | Vt> i| |oo| \w\ || 2 = ^IML 2 - Com- 
bining the estimates for 7 4 , 7 5 and J 6 we get (16.101) . 

Now in order to complete the proof we have to show (16.1 II) . Note that it is useless to multiply 
(16.121) 9 by w since we would obtain a term w 2 div v 2. Thus we adapt again the approach from the 
proof of (|3.2I) and write an expression on a pointwise value of w 2 : 

w 2 (x 1 ,x 2 ) 



f xl f xl Pi 

/ w w s (s, X2)ds = — I —rjrwdivudx 

Jo Jo v\ 

f Xl 1 / \ f xi v\ 2 ^ 

/ — ( w 2 div V2 + wu ■ Vp2j dx — / — t^-w d X2 w dx =: w\ + w\ + 

•J J 1/-^ 



Note that we have p±, v± ~ 1 and thus \/5 > 0: 



/ wjdx <C(\\w\\ L2 \\divu\\ L2 ) <S\\w\\ 2 L2 + C(5)\\u\\ 2 Hl . (6.15) 

J Q 

Next we easily get J^w^dx < E(\\w\\ 2 - j2 + \ \u\ \ 2 H i), and we only have to deal with w 2 . We have 
Jq u>§ dx = Jg 1 [ J p w 2 dx] dx\. Consider the inner integral 

Jp Xl Jp Xi v\ JdP Xl 

The boundary term vanishes and thus 

r (2) 

/ wl < qi^^yiUHlL < JSIMIL- ( 6 - 16 ) 

Choosing for example 5 — \ in (16.151) we get (16.1 II) . what completes the proof of (16.91) . We have 
shown that the solution is unique, and thus comleted the proof of theorem [T] □ 

7 Appendix 

Lemma 9. (Korn inequality) Let V = {v £ H 1 (Q) : (n • v)\r — 0}. Then 3C = C(Q): 

/ 2^T> 2 (u) + vdiv 2 udx>C Q \\u\\ 2 wl . (7.1) 
Jq 2 

The proof can be found in ( ifTOl . Lemma 2. 1) or in ( |f!8)| . Lemma 4). 



Lemma 10. (Helmoltz decomposition) There exists a couple of functions (0, A) G (W 2 ) 2 such 
thatn- V L A\ T = 

u e = V0 + V ± A. (7.2) 

Moreover, 

\w$ + \\A\\w2 <C\\u\\wi. (7.3) 
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The proof can be found in 0. The last auxiliary result we need are the interpolation inequal- 
ities. 

Lemma 11. (interpolation inequalities): 
Ve > 3C(e, p, Q) such that V/ e W^(g): 

||/|| ip <e||V/|U P + C(e,p,g)||/|| i2 . (7.4) 

andVr] > swc/z f&a? ~ + 77 < 1, Ve > 3C(?7, e, p, Q) such that V/ G 

ll/ll w g/^<e||/l|HS + C(c,P,G)ll/llv (7-5) 
Proof. Inequality (17.41 ) results from the following inequality ([1], Theorem 5.8): 

ll/lk<^||/||^li/Hi; e (7.6) 

for each 2 < p < 00, where 6 = "^~ 2 - > and K = K(p,Q). Using Cauchy inequality with e we 
getHH 

The inequality (17.51) for the functions defined on R n is a well-known result from the theory 
of Besov spaces ( GUI ). It can be extended for the square domain due to its symmetry. □. 
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